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Abstract 

We propose an operational method for the solution of differential equations involv- 
ing vector products. The technique we propose is based on the use of the evolution 
operator, defined in such a way that the wealth of techniques developed within the 
context of quantum mechanics can also be exploited for classical problems. We discuss 
the application of the method to the solution of the Lorentz-type equations. 

1 Introduction 

Operational methods provide powerful techniques to solve problems both in classical and 
quantum mechanics. The distinctive feature of these tools is their versatility and the possibil- 
ity of exploiting them in absolutely different contexts, from the time dependent Schrodinger 
problems to the charged beam transport in accelerators. Differential equations have been 
the primary motivation for the introduction of these techniques. Lie introduced the algebras 
bearing his name within the framework of a program aimed to clarify the reasons why only 
restricted families of ordinary differential equations (ODE) can be solved by quadrature. 
Operational methods have become "popular" in applied science for their wide flexibility and 
have stimulated the development of new computer languages, useful for symbolic manipula- 
tion. 

In this paper we go back to the solution of some differential equations involving vector 
products and we will further develop the point of view suggested in Ref. fl], by presenting 
an analysis which includes a variety of problems often encountered in applications, including 
classical time ordering techniques, which are not widespread known as they should. 

We will start our analysis with the following Cauchy problem 

^5 = fx5, SU=o = So, (1) 

almost ubiquitous in physics, from classical mechanics to nuclear magnetic resonance. Al- 
though this equation can always be written in a matrix form, we will develop our consider- 
ations using the vector notation and the properties of the vector product because they are 
more concise and more insightful from the physical point of view. We will assume, for the 
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moment, that the torque vector T is not exphcitly time dependent, so that a straightforward 
apphcation of the evolution operator formahsm yields 

S{t) = U{t)So, U{t)=e'^ (f>(0) = i), (2) 

where tj{t) is the evolution operator defined in terms of the operator T, whose properties will 
be specified below. The series expansion of the exponential provides the following solution: 

n=0 

where T, called vector evolution operator (VOP), satisfies the following identities: 

f^So = So ■■■ f^So = f X (f X ... ifx{f X So)) ■■■) . (4) 

^ V ' 

n— 1 times 

The use of the cyclical properties of the vector product leads to the following closed form 
for (l) (T= |f|) 

S{t) = cos(Tt) So + sin(Tt) {n x So) + [1 - cos(Tt)] {n ■ 5^,) n (n = f/T) . (5) 

This solution has an almost natural geometrical interpretation, which is recognized as a 
Rodrigues rotation (R. r.) [2]. 

The solution of the non-homogeneous version ofeq. Q 

^^S = fxS + N, S\t=o = So, (6) 

is given by 

S{t) = (lit) (^So + ^ dt'U-\t')N^ , (7) 

which is the solution of the first ODE with the evolution operator in place of the integration 
factor. 

The results we have obtained so far are a more concise form for the well known solutions of 
equations of the type (|6]), often encountered in the study of problems involving the Coriolis [3] 
and Lorentz fi] forces. In the next sections we will discuss specific applications of the outlined 
formalism and we will see how it may provide further progress when applied to actual physical 
problems, including time-dependent or spatially non-homogeneous fields. 



2 The Lorentz equation of motion 

The non-relativistic dynamics of a particle with mass m and charge e under the combined 
influence of static electric and magnetic fields, is ruled by the Hamiltonian (c = 1): 

H = ^{p-eA)^ + , (8) 
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where p is the canonical momentum. In the static symmetric gauge, the vector and scalar 
potentials are given by [s] 



A = - B X r , 
2 



$ = -E ■ f . 



(9) 



The equation of motion for the mechanical momentum 

^ ^ e ^ ^ 
71 = mv = p — -B X r 

derived from eq. ([S]), is the Lorentz equation: 

d ^ 



dt 



X V + Q 



(10) 



where we have introduced the following vectors 



n = B 
m 



Q = -E 
m 



This equation has the same form of eqs. ([6j) {T = — ^l, N = Q) and its solution is given by 
eq. ^ with U{t) = e"*^. Therefore, we get (n = B/B): 



v{tj = cos{Uctj Vq + Q + [n ■ ijn — n x m 



UJr. 



where ujc = ^ = eB /m is the cyclotron frequency, and we put 

sin(a;ct) 



i = [1 — cos(a;ct)] Vq + 



Q , m = — ^i{t) 

Ur dt 



A further integration, with respect to the time, yields the position vector, which reads: 

^/ N ^ sm(uct) ^ 1 — cos(uct) ^ 1 ^ 7* 

r{t) = fo + \ ' vo + ' Q + {n-o)n - —nx e , (12) 



OJr. 



where 



sin((X'ct) 



Vo + 



t^ 1 — COs((X'ct) 

2" 



Q. 



(The quantity = vq/uJc is called Larmor radius). The trajectories of an electron for 
different fields configurations are reported in Fig. ([T]). 

It is interesting to stress the contribution n x m, that appears when the initial velocity 
and/or the electric field are not parallel to the magnetic field. It is given by: 



n X m = sm Wet — 2sm 

\ c J ^ 



2 I 



Vd 



ExB 



where Vd is the drift velocity . In physical terms it can be understood as the component 
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y 



y 



Figure 1: Trajectory followed by an electron for different configurations of electric and mag- 
netic fields and initial velocity: (upper panel, left) vq = {vq, 0,0), B = (0, B/\/2, B / \/2),E = 
(upper panel, middle) vq = {vo/V2,0,vo/V2),B = {0,0, B),E = {0,0, -E) (up- 
per panel, right) Vq = {vo/V2,0,Vo/V2), B = {0,0, B),E = {0,0,— E) (lower panel, 
left) vo = {vo/V2,0,vo/V2),B = {0,0, B),E = {0, E/V2, E/^/2) (lower panel, right) 
vo = {V0/V2, 0, vo/V2),B = (0, 0,B),E = (0, E, 0). 



of the velocity which allows the balance between electric and magnetic forces in the Lorentz 
equation of motion. In fact, if we decompose the velocity as v = u + V, and impose that 

eB xV = eE, (14) 

solving for V, we recognize that, in the case B ■ V = 0, this vector coincides with the drift 
velocity vector defined in eq. (13). It is evident that the presence of a further force acting 
on the particle will induce an analogous drift. By replacing E with E + E' in the Lorentz 
equation, we obtain the composed drift velocity 



(^E + E'^ X B 



V, = ^ ^ . (15) 

The role played by the drift velocity can be further stressed by noting that it can be associated 
to the conservation of the so called pseudo-momentum vector [Sl-fT] 



U = m {v-Vd) + e [B X f - tE) , (16) 



which is ensured by eqs. (10), and has been used as an a posteriori check of the correctness 
of our computations. 
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3 Lorentz-type equations and damping 



It is evident that terms of non-electric/magnetic natur^can be added to eq. ([s]). An example 
is represented by the simphfied Drude-hke models [s] , where is introduced a term depending 
on the velocity and a relaxation time and that may be associated with electromagnetic-type 
interactions. From the mathematical point of view the problem to be treated is the search 
of the solution for the following vector differential equation 

^iJ=-nxv + Q--v. (17) 
dt r 

The presence of a velocity-dependent contribution does not modify the procedure described 
before and the solution is given by eq. ([T]) written for the following evolution operator 

U{t) = exp\-t(- + n 



T 



and N = Q. 



The Lorentz equation ( 10 ) is mathematically equivalent to the one describing the falling 
of a body under the influence of the Coriolis force. It has been stressed, within the framework 
of the so called gravito-magnetic theories lo], that effects like the Foucault rotation can be 



associated to a vector potential of the type reported in eq. ^ and the relevant analysis 
has been conducted using the perspective of parallel transport and covariant derivative 10 



The supposed existence of such a potential has given the opportunity of developing further 



speculations as those associated with a possible observation of the Aharonov-Bohm 11 



effect involving the Coriolis vector potential. This phenomenon, suggested by Aharanov and 
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Carmi 12 , has been observed experimentally using the techniques of neutron interferometry 



For the above reasons we will treat in detail the solution of the equation of motion of a 
body falling on the Earth under the action of the Coriolis force. If we include the effect of 
a velocity-dependent force, due, for example, to air friction, we can write this equation as[^ 

—V = — 2uxv — 7]v + g, {tu = un) , (19) 

yJiv 

where u is the angular velocity of Earth rotation. Also in this case, with an obvious redefini- 
tion of the vectors, the formal solution is given by eq. ([T]). The explicit form for the velocity 
vector is: 

{/(t) = e"''* cos{2ujt)vo + {a - e"''* [a cos{2ujt) - b sin(2a;t)]} g 

+ {n-I)n — n X m , (20) 



^Wc mean that the magnetic and/or electric field do not appear explicitly in its expression. As shown in 
the following, genuine non-electric/magnetic effects, like the gravitational term, can also be included. 

^The term associated to the centrifugal force has not been included. We will comment later on the 
possibility of including this type of contribution in the formalism discussed in this paper. 
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where 



e"''* [1 - cos{2ujt)] vo 



— a 



V 



a cos(2u;t) + b sm{2ut) 



m = e sm{2ujt)vo + {6 - e [a sin(2c<;t) + b cos{2ujt)]} g , 



a 



rf + Auj"^ 

while, for the position vector one has 
r{t) 



rf + 4a;^ 



fo + {a - 
+ [at 



e [a cos(2wt) - b sin(2wt)]} -yo 

+ b^ + e^*?* [(a^ - b^) cos(2wt) - 2ab sin(2wt)] } g (21) 



+ {n ■ o)n — n X p , 



(22) 



where 



+ 



P 



a — e 

1 
V 

+ e-^* 



-Tit 



a t 



a cos(2a;t) + b sin(2a;t) 



^^0 



+ a" 



(23) 



b"^) cos(2wt) + 2ab sin(2wt) 



{b - e"''* [a sm{2ut) + b cos{2ujt)]} Vq 

+ [bt - 2ab + e"''* [(a^ - 6^) sin(2a;t) + 2ab cos(2a;t)]} g 



In Fig. ([2]) we show the trajectory described by a falhng body under the influence of Coriolis 
and friction forces, at a latitude of 45° N. The time dependence of the component of velocity 
along the for different values of latitude is shown in Fig. ([s]). 




Figure 2: Trajectory of a body falling under the influence of the Coriolis force and the air 
friction (A = 45° N). The altitude z has been multiplied by 0.01, while the x coordinate has 
been amplified by a factor 500; the body is initially at rest. 
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Figure 3: Time dependence of the x-component of the velocity for different values of the 
latitude. 



We expect that after some time the motion will be dominated by the so called limit 
velocity occurring in any problem characterized by a damping term due to friction, and 
that is reached when the total force acting on the moving body is zero. By imposing this 



condition to the case of eq. (19), we find 



2uj X V* + 7]v* — g 







(24) 



This is a kind of algebraic equation having the velocity vector as the unknown quantity. By 
applying the method developed in this paper, we find for this equation the following formal 
solution 

v* = ^wt9, (25) 
and taking into account the Laplace transform identity 



1 

A 



dse" 



(26) 



valid also if A is an operator, we obtain the following expression for the limit velocity: 



dse-^'^e"'^^ 



772 + 4cj2 



V 



(27) 



This expression yields the generalization of the limit velocity occurring in one dimensional 
motions, but it also generalizes the definition of the drift velocity. In Fig. ^ we have 
reported the modulus of the velocity vector as a function of the modulus of the position 
vector. It is evident that after a very rapid growth, the velocity of the system reaches a limit 
value consistent with eq. (27). In Fig. ([s]) we have reported the components of the velocity 
vector as function of the time. The effect of the Coriolis force is evident from the appearence 
of the components and Vy, which have been significantly amplified to make them visible 
on the same scale for the component Vz- 
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The results obtained in the present and the previous sections are vahd if the vectors 
defining the VOPs do not depend on time and space coordinates. If is not the case, we have 
a completely different phenomenology, that will be discussed in the next sections. 




Figure 4: The modulus of the velocity vector as a function the modulus of the position 
vector; the dashed line refers to the modulus of the limit velocity vector (same parameters 
of Fig. ([2])). The numerical and analytical solutions yields the same result, thus proving the 
correctness of the procedure. 



4 Field inhomogeneities and the Lorentz-type equa- 
tions 

Before getting into the specific details relevant to the solution of the equation of Lorentz- 
type for space and/or time dependent fields, it is important to remind some important issues 
regarding operational ordering problems emerging from the non-commutative nature of the 
vector product. 

We will illustrate the associated difficulties considering, as an example, the Lorentz equa- 
tion in absence of electric field and with a magnetic field consisting of two (non parallel) 
components B = Bi + B2. We may wonder whether the Lorentz equation can be solved in 
such a way that the solution refiects the above decomposition. The evolution operator can 
be expressed in terms of the VOPs VLi^2 associated with the two components of the mag- 
netic field, as f/ = exp \t{VLi + 1^2)], but, as a consequence of the non-commutative character 
of the vector product, the exponential function does not possess the semi-group property 
exp(yl -|- 5) = exp(yl) exp(5) = exp(5)exp(A) [M]. More sophisticated disentanglement 
procedures of the exponential are in order, as, for example, the Zassenhaus identity [Ts] , 
which provides a formula yielding a disentangled expression for the exponential exp (A + B) 
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Figure 5: The components of the velocity vector as a function of time (same parameters of 
Fig. (|4l) 



in terms of successive commutators of A and B, namely 



e^e^e^^e^^ 



(28) 



with 



Zi = - 



1 

2 L 



A, B 



1 
6 



5, 



A, B 



A. 



A, B 



The evolution operator containing the sum of two non parallel VOPs cannot be naively 
disentangled and, accordingly, the associated dynamics cannot be simply expressed as two 



successive R. r.. The identity (28) can be applied and we will illustrate its usefulness by 
keeping only the Zi correction 



f2i , 



and, thus, writing the evolution operator as follow^ 



(29) 



The truncation of the Zassenhaus formula at the first commutator holds for quasi-parallel 
vectors and/or for small times, i.e. when the following inequality is satisfied 

t^\ni X ^l2| < 1 , (30) 

and, therefore, the approximate solution for the velocity vector can be written as 



v'{t) 



(31) 



i exp(|f2i) exp(iii2) exp(|fii) yields an inte- 
gration more accurate than that provided by eq. (291 because it is of the order 0{t^). The inclusion of 



^The symmetric split disentanglement exp[t(f2i + 0,2 



further orders in the Zassenhaus expansion may gives a better approximation. The symmetric split provides 
an easier interpretation in geometrical terms, since it can be understood as three successive R. r. 
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where, on account of the condition (30) and eq. ([s]), v' can be written as 

v\t) ~ e-*'(^i^^2)t;o ~ t/o - Wc,it^c,2t' (^1 X na) x Vo , (32) 

with 

0Jc,k = , Uk = — (fc = 1,2) . (33) 

m Bk 

The successive action of the exponential operators is that of providing two consecutive R. r. 
of the vector v' . 

For example, if we consider the motion of a charged particle under the action of the 
terrestrial magnetic field, gravity and Coriolis force, we should write the equations of motion 
as ^ 

m — = — (e Bt + 2muj) x v + mg , (34) 

where B^ denotes the terrestrial magnetic field. According to the previous discussion, the 
above equation can be solved by introducing a kind of equivalent magnetic field vector given 
by 

-* 777 

B* = Bt + 2 — u , (35) 

e 

thus getting for the associated drift velocity the expression 



m gx B* 

"^=7^B^' ^^^^ 



that give rise to a current flow orthogonal to the gravity force line and to the direction of the 



equivalent magnetic force (35). However, if we are interested to disentangle the magnetic and 
Coriolis components we can follow the just outlined procedure. By assuming that Coriolis 
and Lorentz force vectors are quasi-parallel, i.e. Bt x a; ~ 0, we obtain a first correction 



induced by the combined action of the fields given by (see eq. (32)) 



~ uic,T t sm A sm x , (37) 

where Uc,t is the cyclotron frequency associated to Bt, A is the angle between the vectors 
Bt and cJ, and x the angle between the vector Bt x uj and the initial velocity vo- 

Let us now assume that the magnetic field is not homogeneous, i.e. it exhibits a depen- 
dence on the transverse coordinates. This modifies the dynamics of the charge undergoing 
the Lorentz force effect since during its motion the particle experiences space regions with 
different magnetic field intensity and orientation. Furthermore, by assuming that we can 
choose a (small) finite time integration step 6 during which the fields remains constant, 
the velocity and position vectors can be followed, step by step, by means of the following 
equations (see eqs. (u\ [l2| with Q = 0) 

Vk+i = cos{QkS)vk + [nk-h{S)]nk - Uk x mk{6) , 

ffc+i = fk + ^^"(f^fe _^ . Ok{S)]nk - T^nfc X 4(5) , (38) 
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where the index k corresponds to successive integration steps, and 



£k{S) = [1 - cos{QkS)] Vk Ok{6) 



Qk 



Vk 



Bk 



B{fk) 



m 



We report in Figs. ([6| [7j) an example of motion in a non homogeneous magnetic field, where 
is evident that the main effect of such a coordinate dependence is the appearance of a drift 



velocity contribution 16 . 




Figure 6: Electron motion in a smoothly varying magnetic field lying in the {x,y) plane; 
B = Boix,-y,0),vo = vo{0,l,0). 





Figure 7: Projection of the motion of Fig. ^ in: {z,x) plane (left), {y,z) plane (middle) 
and {x,y) plane (right). The scales are not same (the y-coordinates have been reduced by a 
factor 100) to evidentiate the displacement effect. 



The physical origin of the drift is in the fact that an increase or a reduction of the magnetic 
field implies a corresponding reduction or increase of the Larmor radius, and, therefore, in 
one period the orbit described by the particle is no more closed and the particle, according to 
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the sign of its charge, drifts along the varying field direction. A well known, straightforward, 
calculation allows the evaluation of the drift force under the assumption that the field does 
not vary significantly over a Larmor radius. The magnetic field dependence on the vector 
position yields 



B{fo + 6r) = e'""-^ B{fo) - Bo + (5f ■ V)5o (Bo = B{fo)) 

The extra-contribution to the Lorentz force is 

F = -e < (5f - V)5o XV > , 
where the average is taken over one cyclotron period. In the case Bq 



(39) 



5 ~ i?0! from eqs. (|ll 
F 



12|) one obtains 

eBo 



n X ri) ■ V 



Vt 



Vt = n X vq 



(40) 

Vq, and assuming 
(41) 



and the associated drift velocity is evaluated according to eq. (15). Other types of drift can 
be included, but the procedure remains the same. 



5 Time-dependent fields and Lorentz-type equations 

The electric and magnetic fields in the Lorentz equation of motion may be time-dependent. 
The solution of the problem in the most general cases presents various difficulties associated 
to the fact that the field vectors evolve in time and may not be parallel to themselves at 
different times. This situation is reminescent of what occurs in quantum mechanics when the 
hamiltonian is time-dependent and, thus, does not commute with itself at different times. 
In this case the solution of the problem demands for the use of the ordering methods, which 
can also be exploited for the present problem. 

However, let us start with the case in which the fields are time-dependent but their evolu- 
tion implies only a variation in the modulus but not in the direction. Under this hyphotesis 
time-ordering techniques are not necessary, but the problem deserves some comments re- 
ported below. We will indeed solve the Lorentz equation of motion for a charge moving in 
two mutually orthogonal electric and magnetic fields, with a sinusoidal time dependence, i.e. 
we assum^ 

E = Eq sin((X't + ip) , B = Bq sin^ut + ip) Cy . (42) 



According to eq. (42) the field vectors remain mutually orthogonal and parallel to them- 
selves at any time and, therefore, the solution of the equation is given by eq. ([T]) with the 
evolution operator 

U{t) = e*W (43) 



^With this choice the Lorentz equation we are going to study are relevant for the motion of a charged 
particle in the field of an electromagnetic wave. We have not used the wave-like form sin(fcz — ujt) because 
if we limit the analysis to the non-relativistic case, kz <C ujt (see also ref. [6]). 
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where 

- 1 - e ^ 

= — \cOs(uit + (p) — cos 09] fin ( fio = — Bn Cy 

00 \ m 



Albeit trivial from the mathematical point of view, we report the solution in the simplified 
case E = because it presents some aspects useful for next developments (n = Qq/Qq): 

v{t) = cos(<l>(t))-uo + [1 - cos($(t))] {n-vo)n - sin($(t))n x vq . (44) 

By assuming (f = 7r/2, n ■ vq = 0, and using the Jacobi- Anger expansion [TTj, one obtains 

00 

cos(xsin6') = Jo{x) + 2 cos(2?26') J2„(x) , 

n=l 

00 

sin(xsin^) = 2 ^ sin[(2n + 1)6] J2„+i(x) , (45) 



n=0 



we obtain, from eq. (44), the following expression for the velocity {( = ^0/^) 

v(t) = Jq{C)^o " 2Ji(^) sm{ut)nxvo 

00 

+ 2 ^ {cos{2nujt) J2n(C) ^0 - sin[(2?2 + l)ut] J2n+i(C) nxvo} , (46) 



n=l 



which shows the interplay between the cyclotron frequency and the frequency of the oscil- 
lating magnetic field. 

The inclusion of the electric field implies some additional computational problems. Again 
in the case ip = 7i/2, we write the inhomogeneous term of the solution as follows 



w 



(t) = e*W rdt'e-*(*')g(t') 
Jo 



g-sinM)C ^ ^^/gSi„M')C cos(a;t')<5o ( <5o = ^ 4 ) • (47) 



This integral can be treated in different ways, and, to have an idea of the mathematical 
problem one may face with, we choose to carry out the integration using the Bessel function 
expansion of eq. (45). We get: 

wit) = - 2^ 2^ exp U^— F„(t) Jfc,„ (48) 

fe=— 00 n=— 00 ^ ^ 

where 
and 

\n = Jk{^C)Jn{-^C)Qo■ (49) 
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The last vector can be specified either in terms of the series expansion for the Bessel func- 
tion^ or by the use of the integral representation 17 



that allows to write 

— * 

Jk.n 



1 

4^ 



1 

2^ 



de 

— 7T J — TT 



dxe 



-i(fce+nx) -{sin0-sinx) C 



(50) 



(51) 



We have mentioned this specific problem to give a very first idea of the problems associ- 



ated with the solution of equations of the type (10), where the torque and inhomogeneous 



vector are explicitly dependent on the integration variation, but also because it has inter- 
esting implications for the understanding of the role played by the Poynting vector in the 
dynamics of charged particles moving under the combined action of mutually time dependent 
orthogonal fields. The method of solution we have proposed, having an intrinsic vector na- 
ture, can be ideally suited to treat the question addressed in and clarify the link between 
the Poynting vector and the drift velocity term given in eq. (13). This aspect of the problem 
will be treated elsewhere. 



5.1 Operational Methods and time ordering techniques 

Before considering problems requiring time ordered products, we will present a method 



of solution involving an elaboration of the Heaviside operational method 18 . To better 
appreciate the usefulness of the procedure we consider the case in which the particle is 
initially at rest and only the electric field is a function of time while the magnetic field is 



^(^) = ( :^ + ^ ) <5(t) , (52) 



static. We can write the formal solution of eq. (10) in the following way 

-1 

and the use of the Laplace transform methods leads to 

v{t) = ds exp |-s (^A + I Q{t) . (53) 

In this integral the exponential can be straightforwardly disentangled because the VOP and 
derivative operators appearing in its argument commute between them. Therefore, from eq. 



(11), under the further assumption that the electric field differs from zero only for t > 0, one 



has {n = VL/VL) 

v{t) = I dt'e^''-'^^Q{t') = c + {n- f)n + fix s, (54) 







Jn{z) — J2k=0 fc!(n+fe)! (2/2)"+^*^ 
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with 



c = [ dt' cos [n{t' - t)] Q{t') , s = f dt' sin [n{t' - t)] Q{t') , 

Jo Jo 

f = [ dt' {1 - COS [Q{t' - t)]} Q{t') . (55) 
Jo 

The case in which the vector VL varies with time, not only in modulus but also in direction, 
in such a way that 

n(ti) X Q{t2) ^ , (56) 

implies that the VOPs associated at different times do not commute. The solution of our 
problem cannot be obtained using a straightforward integration of the time dependent part, 
i.e. it will not be sufficient replace in the Rodrigues rotation tCl with drf2(r). 



From the geometrical point of view the condition (56) states that the torque vector 
is no more parallel to itself at different times, or that the corresponding matrix equation 
is expressed in terms of an explicitly time-dependent matrix, not commuting with itself 
at different times. The situation is clearly reminiscent of what is occurring in quantum 
mechanics where the solution of Schrodinger problems requires a time ordered expansion 



of the evolution operator, like the Dyson 19 or Magnus 20 expansion. We will treat the 



problem by exploiting the theory of path ordered exponential |2T|. 

From the mathematical point of view the path ordered exponential function is defined 
in non commutative fields and is equivalent to the exponential function in a commutative 
field. We define, therefore, the following ordered exponential with respect to the ordering 
parameter t: 

oe[f]{t) = (exp\ f dt'fit')]] (57) 



where the symbol ()+ denotes the Dyson time ordering operator for the element T{t) of an 
algebra with a non commutative product o. The ordered exponential can be defined in many 
different ways. Here we use the differential equation 

^oe[f](t) = f{t) o oe[f](t) , (oe[f](0) = 1) . (58) 



The solution of eq. (58) can be written in terms of the following series 

oo 

oe[f](t) = i + ^ oe„[f](t) (59) 



n=l 

with 

ft 



oe^[f]{t) = - I dt'f{t') o oe„_i[f](t') , oeo[f](t) = 1 . 
n\ Jo 

The above solution is clearly recognized as an ordinary Dyson expansion ^^I9j, which is not 
fully satisfactory because it is a perturbative series which does not ensure properties, e.g. 
the conservation of the norm of a vector, holding also for time-dependent vectors. Different 



15 



expansions, preserving at any order the norm, can be employed as, for example, the already 
quoted Magnus expansion (see also ref. |2l]), which can be written as follows 



oe[T]{t) = exp<j j dt'T{t') 



1 



t rti 
O , / dt2 

-h Jo 



(60) 



where the dots refer to higher order commutators, not reported here. Let us note that in 
the case of Lorentz equation the non commutative product is the vector product (o = x), 
and, therefore 



f{ti),f{t2) = 2T(ti) o T(t2) 



(61) 



According to this result, the solution of the Lorentz equation (10) with f2 and Q depending 
on time is given by eq. () with the evolution operator 



U{t) = e-^W+^W 



where the following notation has been introduced 



drfi(r) 



A(t) 



t ph 

dti / dta 
Jo 



Q{ti), Q{t2 



+ 



(62) 



(63) 



As for the correction A(t), if we assume that the modulus of the vector Q remains 
constant, one has 



2 sin 9i2 u , 



(64) 



where 612 is the angle formed by the two vectors and u is the versor pointing in the direction 
orthogonal to the plane defined by fl{ti) and ^(^2)- For sufficiently small time differences 
and for adiabatic changes, we expect that ^^12 = a;(t2 — ^i) ^ 1, and therefore 



A(t) ~ -Q^ut^u 
6 



(65) 



Higher orders corrections can also be included but calculations becomes more and more 
cumbersome. 

The previous discussion has been developed on purely mathematical grounds. As an 
example of application, we can consider the motion of a particle under the influence of a 
magnetic field with a slowly varying component along the ^-axis and a constant y-component, 
i.e. the magnetic vector changes its direction and it is not parallel to itself at any time. The 
lowest order corrections in the Magnus expansion allow the inclusion of the effect of the 
slow field evolution. If we assume that the variation is adiabatic over one gyration period. 



the correction can be evaluated by means of eq. (65), with a frequency u, assumed to be 
constant during this time, given by 



ByBz 



(66) 
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6 Second order Lorentz equation 



Before entering the specific topic of tliis section we will discuss the application of the oper- 
ational method developed in the previous sections to the solution of the following evolution 
equation 

^-5 = f X 5 + AT X ff X 5) , ^|t=o = 5^o, (67) 



dt 

which is an evolution-type vector equation, with a further contribution associated with a 
double vector product of the torque vector. The associated evolution operator writes 

U(t) = e*^ + ^*^' (68) 

which involves linear and quadratic VOPs. It resembles the generating function of two 
variable Hermite polynomials i?„(a, b) [23] and can be therefore expanded in series according 
to the identity 

'^'^''^ = E f^^n(a,6) H^{a,b) = n\ 7^^^^ • (69) 

The use of the operator T as expansion parameter yields the following series for the evolution 
operator 

00 _ 

ij{t) = J2 -,H4t,Xt)f-. (70) 
^ — ^ n! 

n=0 

Even though slightly more complicated than an ordinary exponential expansion, we can 
again take advantage from the cyclic properties of the vector product to reduce it to a kind 
of Rodrigues rotation, by defining the following cos- and sin-like functions 

n=0 n=0 

The form of the solution is therefore exactly that given in eq. ([S]) with the replacements 
cos — )■ Ch, sin — )■ Sh. In the case in which the evolution operator 

Uit) = expjj] Afctf'^l (72) 

we can exploit the same procedure involving higher orders Hermite polynomials, as we will 
discuss in the concluding section. 

Let us now discuss the inclusion of radiation correction effects, which are usually incor- 
porated in the classical Lorentz equation by means of a second order time derivative term, 
according to the following expression [2] 
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where 



dt 



v\t=o 



2 ro 

3 c 



In this case equation we have the contribution of an extra term that acts anti- 
damping giving rise to the so called runaway solutions. The physical content of this equa- 
tion is well known and will not be commented here. We will limit our analysis to its math- 
ematical aspects, which have some elements of interest, since this equation represents a 
non-homogeneous second order differential vector equation and we can use an extension of 
the previously outlined method to write the relevant solution. For future convenience we 
factorize the operator acting on the velocity vector as follows 



d 
dt 



A, 



dt 



- A_ 



1 



Q 



(74) 



where 



A 



± 



1 

27 



:i ± a) 



a 



1 + M 



In the hypothesis that Vq = do = 0, we can write the formal solution of eq. (73) using a 
generalization of eq. (52), namely 

m = -z d-A^] ' Q{t) 



d 
dt 

1 



r(A+ - A_ 




A, 



-1 



dt 



-1 



Q{t) 



from which, by using the Laplace transform identity (26) and assuming that Q 
we obtain 

ft 



(75) 



for t < 0, 



v{t) = - 



a 



dee-«^+g(e) - e 



tA. 



dee-«^-g(0 



(76) 



We note that 



1 

a 



00 -s{1 + 4tQ) 

ds — 







and, using the Newton series expansion for the operator A±, 
operators: 



(77) 

we get for the exponential 



exp 



k=0 



1/2 
k 



(7J 



which can be written in terms of a Rodrigues rotation involving the previously quoted sin- 
and cos-like functions expressed in terms of the higher orders Hermite polynomials (see eq. 
(0)- 

According to the previously outlined steps we know how to handle the formal expression 
given in eq. (78) to get an explicit solution for our problem. Let us now consider only the 

hose 

m 



homogeneous part of eq. (73), whose formal solution reads 



ci + e*^- C2 



(79) 
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where Ci,2 are constant vectors linked to the initial vectors by the relations 



(i+ - i_)ci 

(i+ - i_)c2 



-A_ Vq + So 

A+ vo - ao . 



(80) 



The action of the exponential operators on the constant vectors can be defined according to 
the previous prescriptions. 



Since eq. (73) is a second order equation, we can cast it in the form of a matrix equation 
as follows 

^ Z = MZ + K , (81) 



dt 



where 



M 



1 

r 



r 
tl 1 



K = 



T 





Q 



12) 



The evolution operator e**^ associated to this equation is given by the exponential of a 
2x2 matrix. Since the matrix M does not depend on time, standard means, e.g. the 



Cayley-Hamilton theorem 22 , can be used to cast it in the form reported below (cr = t/2r) 

uit) 



tM 



4 sinh(crd) + cosh(cra) 2^ sinh((jd) 

2 ? sinh((T a) 4 sinh(cr a) + cosh((T a) 



and, thus (Zq 



\t=o) 



Z = U{t)Zo + / dt'U-\t')K, 
Jo 

that is the same solution as before, written in matrix notation. 



83) 



The structure of the previous formalism may rises some confusion. As an example, let 
us consider the definition of the norm of Z: 



{v a) 



^5) 



where the dot represents an ordinary scalar product. It is evident that, according to such 
a definition, the norm is not preserved. This is due to the rotation matrix, which is not 
norm-preserving, and to the runaway mechanism, associated with the anti-damping term e'^. 
Such a term is, from the mathematical point of view, not particularly significant and can 



always be eliminated by means of a Liouville transformation 17 



Let us consider again the following second order vector equation 

d^ . d 

— s + A — s + Bs = Q 
dt"' dt 



(86) 



where A and B are VOPs. An equation of this type is met in physics in the study of Coriolis 
problems, including the effects of the centrifugal forces, and presents an extra difficulty 
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associated with the presence of the vector operator in the damping term. Eq. (86) can be 
reduced to a Liouville standard form by setting 



s{t) = exp --i u{t) 



(87) 



and, if [A, B] = 0, we can write the equation for the vector u as follows 



dt2 



u+ [B 



4 



A' \ u = 0. 



The transformation provided by eq. (87) is geometrically interpreted as a R. r. of the vector 



u, induced by the torque vector associated with the VOP A. The solution of eq. (88) can 
be obtained using the procedure illustrated before. 

The previous results (see also the first paper in Ref. |l]) can be exploited to develop a 
numerical code for the motion of bodies under the influence of gravity, Coriolis and centrifugal 
forces. 



7 Relativistic effects 

In the previous sections we have discussed the motion of charged particles in electric and 
magnetic field without considering any relativistic correction. This can be easily accounted 
for by rewriting the Lorentz equation in the form 

mo ^ i-fv) = e [-B x v + (89) 

where mo is the electron rest mass and 7 is the relativistic factor. The integration of this 
equation appears problematic, even for constant and homogeneous fields, since E induces a 
non conservation of the modulus of the velocity. The relativistic factor 7 is no more constant 



and a further equation, specifying its time dependence, should be coupled to eq. (89). By 



introducing the vector A = ^yv, eq. (89) can be rewritten as (with the usual means for the 
vectors Q and Q) 

= --nx A + Q , (90) 
at 7 

By multiplying both sides of this equation by the velocity, we obtain 

A-^A = A-g. (91) 

and, thus (Aq = A(t = 0)) 

A\t) - Al = 2 [ dt'Q ■ A(t') , (92) 
Jo 

that provides the relativistic kinetic energy variation due to the interaction of the charge with 
the electric field. As already stressed, the presence of the relativistic factor prevents us from 
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the possibility of finding an analytical solution using the tools developed so far. However, a 
fairly straightforward integration scheme can be used adopting the same iterative procedure 
described in sec. [i] (see eqs. (38)). In this case the solution is obtained at each time 
step in terms of eqs. (11), in which Q is replaced by f2/7„_i. An example of solution of the 
relativistic problem is shown in Fig. (^ where we have plotted the modulus of the velocity as 



function of time using the solution given in eqs. ( 11 ) and including the relativistic correction. 
It is evident that for low initial velocities the two solutions are undistinguishable, while the 





Figure 8: Modulus of the velocity as function of time: comparison between the non relativis- 
tic (left) and the relativistic (right) case for vq = fo(l, 0,0); B = Bo{0, 0,1); E = Eo{0, —1, 0). 

differences becomes significant with increasing values of Vq. The oscillations in the velocity 
are due to the combined effect of electric and magnetic fields, but in the relativistic case the 
oscillation period becomes larger, in agreement with the fact that the cyclotron frequency 
decreases with increasing 7. 

We do not include further examples because the goal of this section was just show- 
ing that the analytical solutions could straightforwardly be implemented into a numerical 
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scheme yielding the relativistic treatment. These effects will be thoroughly discussed in a 
forthcoming paper. 



8 Bremsstrahlung effects 

The radiation emitted by an accelerated charge can be evaluated from the Lienard-Wiechert 
integral j4j which yields the energy radiated per unit solid angle and unit frequency as 



dVt dw A'n'^c 



(93) 



with S modulus of the vector 



T 







S= / dt gx (g X /3) exp<'?a; ( t- — ) J> , (94) 



c 



where q denotes the unit vector of the direction along which the emitted radiation is observed 
(see Fig. (|9|), t7 = /3c is the velocity of the charge, and T is the time during which the charge 
effectively experiences the acceleration due to the fields. By assuming that the motion occurs 




Figure 9: Reference frame definition for eq. (94). The magnetic field is directed along the 
z-axis. 



in absence of electric field and under the influence of a constant magnetic field, by using eq. 



exp iiuo [t 



q ■ r 



( [12| (with (fo = 0) it is easy to show that: 

= expjz— ci| exp{«a;c2t} 

exp |z— [cs sin(fit) — ci cos(i7t)]| 



(95) 
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where {n = B/B, f3o = vq/c) 

ci = q - {fix (3o) , C2 = 1 - (ri ■ /3o) (g ■ n) , C3 = {n ■ /3o) (g ■ n) - q- /3o 
and, according to the Jacobi-Anger expansion 



expjz^ [C3 sm{nt) - cs cos(r]t)]} = ^ e'^^' B^ 

k=—oo 

00 00 

fc=— 00 m=— 00 

Moreover, it turns out 

qx{qxP) = ai + be cos{Qt) + 6^ sin(fit) , 
where we introduced the following vectors 

a = in-Po) [{q-n)q - n] , 

be = {q- Po)q - - a, 
hs = q-{nx /3o) - (n x f3o) . 



(96) 



(97) 



(98) 



Inserting these results in eq. (93), one obtains (for notation semplicity we have omitted to 
indicate the argument of the i?-functions) 

T 



S 



where 



expjz^cij J2 + bBr-i + b*Br+i] exp<jz^T|> sinc(^T) , (99) 



h = he + ihs 



r f2 + w C2 . 



According to the previous equation, the spectrum of S presents a series of harmonics with 
frequency 

rVt 

UJr = — . (100) 

1- {q-n) {n- Po) 
In the case of a relativistic electron, we have 

e B Ue,o 



n 



7 mo 



7 



(101) 



^The functions can be written in terms of the cylindrical Bessel functions by exploiting the identity 
In{ix) = J„(x) and the Graf addition theorem 17 , which yields 

%/2 



B, 



\x-iy J 



23 



and we can treat the solenoid as an undulator with period 

2tic 



Xr 



7\o 



7 



If n = q and the initial motion is mainly in the same direction of the field, we find 

Xr 



2 r 7^ 



(102) 



(103) 



If we limit ourselves to the non relativistic case, the inclusion of the electric field does 
not imply any particular computational problem for the radiation integral. A significant 
difference is associated with the lineshape which is not simply the square of the sinc-function 
appearing in eq. (99), but is now given by the square modulus of the function 

at exp \ (prt — u 1 



(104) 



Fig. (10) shows the shift and broadening of the spectral line of the radiation emitted by 
electrons moving in the field of a solenoid. 




Figure 10: Comparison of the emitted radiation lineshape between the case with only the 
solenoid (a) and the case where also an electric field is present (b). 



9 Quantum mechanical aspects and concluding remarks 

The final topic we will treat is the extension of these methods to quantum mechanics. To 
this aim we consider again the Hamiltonian given in eq. (^, but in absence of any electric 
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field. If we assume that tfie magnetic field is directed along the z-axis, we can write the 
explicit form the Hamiltonian as follow^ 

H = — i-ihV - elV 
2m V / 

= {-n^V^ + e^B^{x^ + y'^) - ihe B {xdy - yd^)] , (105) 

^ III 

which can be more conveniently cast in the form (see also [H] and references therein) 



where 



(107) 



The operators Fi 2 satisfies the commutation rule 



(108) 



and (T = f if 2 



-AtT + 2 



\k+l 



2itl (A; = 1,2) 



(109) 



i.e. they exhibit the commutator properties of SU(1,1). The relevant Heisenberg equations 
write in the form of a vector torque equation of the same type of eq. ([T| and the time- 
dependent solution of the associated Schrodinger problem can be obtained using the standard 



Wei-Norman ordering procedures 24 



A different way of treating the quantum evolution problem consists in introducing a 
suitable transform as follows (C = (Ci^ C2, Ca)) 



U{t) = exp 



2mh 



(27r)3/2 J 



dC e ^ exp <( 2 z 



(110) 



tt 



2mh 



(■ {ihV + eA) 



where, from eq. ^ 



{thV + eA)i = thd^ - -By 

{thV + eA)2 = ihdy + ^Bx 
{ihV + 6^4)3 = ihd^ . 



(Ill) 



^Note that in the static symmetric gauge (cfr eq. J9|), V • ^ = 
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In eq. (110), by ordering the exponential in the integrand, for the wave function at later 
time we obtain 



^(f,t) = f/(t)^(r) 



dC expi -C^ + c<;cKiC2 - i \Hrr (^c{y Ci -xC,2) 



, ^ i2fit , ^ 

^{f-J C • 112 

V m 

This relation is a generalization of the Gauss- Weierstrass transform and the study of its 
consequences will be discussed elsewhere. 

It is worth stressing that the case where the direction of the magnetic field is arbitrary 
implies only a slightly more complicated integral transform. Even the inclusion of the electric 
field is not particularly tricky. Assuming that the electric field lies in the (x, ?/)-plane [E = 
{Ex, Ey,0)), we can recast the Hamiltonian operator in the form 

H = — {-ihV - eA] - eE-f 
2m \ / 

e B 

- ^ ^ ^ [(^ ~ '^y - (y - *9x] 

1 zu 
- -mw'^ {al + al) + « ^ y {otx + ay dy) , (113) 



where 



AmE^y eB Uc 

„ = —r^ zu - - 



^•^ eB^ '"' 2 m 2 

This Hamiltonian is that of a multidimensional harmonic oscillator with shifted coordinates. 
The constant term ^mzu'^ (a^. + a^) is just the vacuum field energy redefinition associated 
with the transformation used to move to the shifted coordinate representation. The term 
proportional to [a^ + ay dy) is the quantum counterpart of the drift motion. 

In the paper we have stressed the analogy between Coriolis and Lorentz forces. This 
is more than a formal analogy and the previous considerations can be extended to the so 



called quantum Coriolis states 25 , that are similar to the Landau quantum states 26 
entering in the analysis of the motion of a quantum electron in a classical magnetic field. 
Some aspects of the problem and the possibility of studying them within the context of the 
present formalism, will be the argument of a forthcoming investigation. 

The methods we have developed are flexible, fairly simple and easily amenable for nu- 
merical computation. Their use can also extended to non linear equations like the Landau- 
Lifshitz-Gilbert equation describing the precessional motion of the vector magnetization M 



in solids 27 . Without entering in details, we remind that such equation, in our notation, 
writes 

dtM = -{a + l3Mx){M X H) . (114) 



26 



The quadratic non-linearity creates noticeable difficulties and the equation can be viewed 
as a kind of Riccati vector equation. Assuming however that the conditions of the problem 
allow the definition of a time step in which the variations of the vector M are not large, we 
can use the following solution scheme 

dtMn = PxM^, P = (a + (3 Mr^.ix) H , (115) 

— # — » 

where the vector P, containing the vector M at the previous integration step, is treated as a 
constant torque with the inclusion of an anti-damping term. The solution is essentially that 



provided in eq. (20), and the evolution should be followed by successive steps. 

The points in these concluding remarks have briefly treated just to show the flexibility 
of the method we have proposed and the number of topics which it allows to treat. Most of 
them deserves a deeper analysis, that we will develop in a forthcoming investigation. 
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